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ABSTRACT 

In this paper we survey the theory of wind accretion in high mass X-ray binaries hosting a magnetic 
neutron star and a supergiant companion. We concentrate on the different types of interaction between 
the inflowing wind matter and the neutron star magnetosphere that are relevant when accretion of 
matter onto the neutron star surface is largely inhibited; these include the inhibition through the 
^v^j , centrifugal and magnetic barriers. Expanding on earlier work, we calculate the expected luminosity 

for each regime and derive the conditions under which transition from one regime to another can 
take place. We show that very large luminosity swings (~10 4 or more on time scales as short as 
hours) can result from transitions across different regimes. The activity displayed by supergiant fast 
X-ray transients, a recently discovered class of high mass X-ray binaries in our galaxy, has often been 
interpreted in terms of direct accretion onto a neutron star immersed in an extremely clumpy stellar 
wind. We show here that the transitions across the magnetic and/or centrifugal barriers can explain 
the variability properties of these sources as a results of relatively modest variations in the stellar 
wind velocity and/or density. According to this interpretation we expect that supergiant fast X-ray 
transients which display very large luminosity swings and host a slowly spinning neutron star are 
characterized by magnetar-like fields, irrespective of whether the magnetic or the centrifugal barrier 
applies. Supergiant fast X-ray transients might thus provide a new opportunity to detect and study 
^ . magnetars in binary systems. 

^ ' Subject headings: accretion, accretion disks — stars: neutron — supergiant — X-rays: binaries - 

X-rays: stars 
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1. INTRODUCTION 

GN ! High mass X-ray binaries (HMXBs) consist of a collapsed object, usually a magnetic neutron star (NS), that accretes 
matter from an OB companion star. Mass transfer tak es place because of the intense st ellar wind from the OB star, 
part of which is captured by the collapsed object (e.g. IVerbunt fc van den HeuvellflQQfih . Only in some short orbital 
period system s, the early type star, often a s upergiant, fills its Roche lobe and leads to mass transfer through Roche 
lO ' lobe overflow (|Tauris fc van den~H cuvcl 2006) . Persistent HMXBs accrete all the time and in most cases display X-ray 
CD luminosities in the 10 35 -10 38 erg s -1 range. Many HMXBs are transient systems that remain at low X-ray luminosity 
levels (10 32 -10 33 erg s" 1 ) most of the time and undergo outbursts lasting from weeks to months. During these outbursts 
they display nearly identical properties to those of persistent HMXBs. Transient systems usually comprise a Be star 
donor and relatively long, moderately eccentric orbits, such that the star sits deep in its Roche lobe and stellar wind 
capture is the only mechanism through which mass transfer takes place. The occurrence of the outbursts is likely 
associated to variations in the stellar wind of the Be star, su ch as shell ejection episodes, or build up of matter around 
\ the resonant orbits in the slow equatorial wind component (|van den Heuvel fc Rappap ort 2005). However, there are 
characteristics of the outbursts that are difficult to interpret if accretion onto the neutron star surface takes place 
unimpeded also in quiescence; these are (a) the large outburst to quiescence X-ray luminosity swing (factor of ~10 3 or 
larger) and (b) the presence in a given source of low-luminosity (Type I) outbursts recurring close to periastron and, 
at different times, of high-luminosity (Type II) outbursts that last for several orbital cycles and display little (if any) 
X-ray flux variations associated to the orbital phase. These characteristics of Be transients can be explained if the 
accretion rate (and thus X-ray luminosity) variations that are produced by the stellar wind alone, could be amplified 
by some "gating" mechanism. Since most Be star HMXB transients contain relatively fast spinning X-ray pulsars, 
such mechanis m has been identifi ed with the centrifugal barrier that results from the rotation of the neutron star 
mag netosphere ()Stella et al.lll986l ). 

About 10 transient systems have been recently discovered, which display sporadic outbursts lasting from minutes to 
hours (i.e. much shorter than Be star transients') and reach peak luminosities of ~10 36 -10 37 erg s -1 . These systems 
spend long time intervals in quiescence, with X-ray luminosities down to ^10~ 5 times lower than those in outburst; in 
spite of their association with OB supergiant companions, their behaviour is thus at variance with other persistent and 
transient HMXBs. They define a new class of HMXBs, collectively termed supergiant fast X-ray transients, SFXTs. 
An overview of the properties of SFXTs is given in § [2] 
If accretion onto the collapsed object of SFXTs takes place both in quiescence and outburst, then the corresponding 
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Fig. 1. — Schematic view of a magnetized NS interacting with the inflowing matter from its supergiant companion. All the regimes of 
§[3] are shown, together with the relative position of the magnetospheric radius (solid line), the corotation radius (dashed line), and the 
accretion radius (dotted line). A wavy solid line is used when the magnetospheric boundary at Rm is Kelvin-Helmholtz unstable. In the 
supersonic and subsonic propeller regime convective motions at the base of the atmosphere are represented with small eddies. 



X-ray luminosity swing, typically a factor of ^10 4 -10 5 , would require wind inhomogeneities with a very large density 
and/or velocity contrast (according to the standard wind accretion, the mass capture rate onto the NS scales like 
M w v~ 4 , with M w the mass loss rate and v~ 4 , the wind velocity of the supergiant star. iDavidson fc Ostrikerlll979fl . 
Several authors (|In't Zandl 2005: Le vder et al.1 120071 IWalter fc Zurita Herasl I2007T) suggested the presence of dense 
clumps in the wind of the OB companions in order to attain the luminosity variations of SFXTs. While some 
observations provide evidence for a clumpy wind, the characteristics of such inhomogeneities are still poorly known. 
Numerical simulations suggest that clumps m ay originate from small scale perturbations in the radiation-driven wind 
( Dessart fe Owockl 120031: iPrinia et al.ll2005f). Models involving accretion of clumps are still being actively pursued 
( Negueruela et al.l 120081 : IWalter fc Zurita H cras 2007). The requirement on the density and/or velocity contrasts in 
the wind can be eased if there is a barrier that remains closed during quiescence, halting most of the accretion flow, 
and opens up in outbursts, leading to direct accretion (this is similar to the case of Be star transients). This paper 
develops and discusses gated accretion models for SFXTs. 

In §[3] after reviewing the theory of wind accretion in HMXBs, we describe the different regimes of a rotating magnetic 
neutron star immersed the stellar wind from its companion. In §[3] we discuss the conditions under which transitions 
across regimes can take place in response to variations in the wind parameters. A s compared to pr evious works 
addressing gating mechanisms in transient accreting magnetic neutron stars (notably iStella et al.|[l986f ). we present 
here a more comprehensive treatment of the different physical processes that have been discussed in this context. 

In § [5l we present an application to the different states of two SFXTs and discuss in turn the possibility that the 
outbursts are driven by a cen trifugal or a magnetic barrier. The latter mechanism involves a magnetosphere extending 
beyond the accretion radius (|Bondil ll952). which prevents most of the inflowing matter from accreting; it is discussed 
here for the first time in relation to the activity of transient binary sources. Unless the stellar wind were extremely 
clumpy (as envisaged in other SFXT models), the onset of a barrier inhibiting direct accretion would be required to 
explain the activity of SFXTs. In this case we conclude that if SFXTs host slowly rotating neutron stars (spin periods 
of several hundreds to thousands seconds), then they must possess magnetar-like fields (~10 14 -10 15 G), independent 
of whether the centrifugal or magnetic barrier operates. We summarize our discussion and conclusions in §[6] and § [Tj 

2. THE OBSERVED PROPERTIES OF SFXTS 

SFXTs are observed to exhibit sporadi c outbursts, lasting from minutes to hours, w i th peak X-ray luminosities 
between ~10 36 and IP 37 erg s" 1 (see e.g . iGonzalez-Riestra et al.ll2004t ISidoli et alj|2005t iGrebenev fc Sunvaeyll2005l : 
iLutovinov et al.ll2005t ISguera et ai1l2005l : iMasetti et al.ll2006HSguera et al.ll200fl iGotz et al.ll2007l : ISguera et al.ll200l . 
No firm orbital pe riod measurement has been o btained yet 4 . A recent list of confirmed (~5) and candidate (~6) 
SFXTs is given by IWalter fc Zurita Herasl (|2007f ) . 

Between outbursts, SFXTs remain in quiescence with luminosities in the ran ge ~10 31 -10 33 erg s _1 
(jGonzalez-Riestra et al.ll2004l : lln't Zandll2005l : ISmith et al.l 2006: K ennea fc Campanall200ll . In some cases, very high 
peak-to-quiescence X-ray luminosity swings (factor of ~T0 4 -10 5 ) were seen on timescales com parable to the outburst 
durat ion. Some SFXTs showed also flare- like activity at intermediate luminosity levels (e.g., Gouzalc z-Riestra et al.l 
2004). In the case of IGR JI7544-2619, two states of intermediate luminosity were observed: one before the onset of 
the outburst and the other immediatel y after, with X- ray luminosities ~3 and ~1 decades below the value reached at 
the peak of the outburst, respectively (|In't Z and 2005j). 

The sporadic character of SFXT outbursts, as observed with INTEGRAL, suggested that the duty cycle of these 

4 Only IGR J11215-5952 and AX J1749. 1-2733 showed recurrent flaring activity, with periodic ity of ~165 d and ~185 d, respectively. 
These are interpreted as outbu rsts fr om two systems with unusually long orbital periods (>100 d, Grcbcncv & Sunyacv 2005; Sidoli et al. 
120071: IZurTta Heras et 111120071 1. Thus IWalter fc Zurita Herasl l(2007fi excluded these sources from their SFXT list. 
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sources (the fraction of time spent in a high luminosity state) is small (~0.02-0.002. IWalter & Zurita Herasll2007fh 
However, recent observations carried out with the very sensitive X-ray telescopes on board XMM-Newton and Chandra 
revealed that some SFXTs display flares around a luminosity of >10 34 -10 35 erg s" 1 (i.e. well below the INTEGRA L 
limiting sensitivity) for a large fraction of the time (|Gonzalez-Riestra et al.|[2004t lln't Zand 200]| iTomsick et al.l l2006') . 
Therefore, the indication is that the active phase of SFXT sources (as opposed to true quiescence) lasts longer than 
previously thought, and the duty cycles are of order ~0.1 or higher. 

Optical identifications of S FXTs show that these sources are associated to OB supergiant companion stars (see e.g., 
IWalter fc Zurita Herasl[2007l and reference therein). The SFXT OB companions have typically mass of M„~30Mq, lu- 
minosity of log (L„/L Q )~5-6, mass loss rate of M W =1CP 7 -10~ 5 M Q yr _1 , and wind velocity of 1000-2000 km s _1 . Note, 
that isolated OB stars with log(L < ./L f 7))^5-6 are persistent soft X-ray sources with luminosity around ^-TO 32 erg s _1 
(|Cassinelli et alj|1981t iBerghoefer et al.ll2001h . 

It is widely believed that SFXTs contain sporadically accreting neutron stars. Only little is known about their spin 
period. A coherent periodicity was detecte d at 4.7 s and 228 s in AX J1841. 0-0536 and IGR J16465-4507, respectively 
(|Bamba et al.ll20Mt iLutovinov et al.l 2005). However, IGR J16465-4507 showed only a factor of ^100 luminosity swing 
between quiescence and outburst, so it is unclear whether the sou rce should be considered a transient (in fact it is 
classified as an "inter mediate system" . [Walter fc Zurita Her as 2007). The nature of the companion star in AX J1841.0- 



0536 is still debated ([Halpern et al.ll2004UNespoli et al.ll2007 ). Therefore, these two sources m ight not belong to the 



SFXT class. On th e contrary, in the prototypical SFXTs XTE J1739-302 ([Sguera et alj [20061 and IGR J16479-4514 
([Walter et al.l [2006h some evidence has been reported for periodicities in the ~1000-2000 s range. We assume in the 
following that SFXTs host a rotating magnetic neutron star. 

3. STELLAR WIND ACCRETION 

We investigate here the conditions under which a magnetized neutron star can accrete matter from the wind 
of a massive companion. In the theory of wind accretion in HMXBs, the following radii are defined (see e.g., 
Illlarionov fc Sunvaevi I1975I : IStella et al.lll986f ): 

• The accretio n radius, R a is the distance at which the inflowing matter is gravitationally focused toward the NS 
([Bondilll952D . It is usually expressed as 

i? a = 2GM NS /t4 = 3.7 x 10 1C V cm, (1) 

where vs is t he wind velocity in units of 1000 km s _1 and we assumed that the orbital velocity of the star 
is negligible ([Frank et al.l [2002) . Throughout the paper we fix the NS radius and mass at Rns=10 6 cm and 
Mns=1-4 Mq, respectiv ely. The fraction M ra pt,/M w of the stellar wind mass loss rate (M w ) captured by the NS 
depends on R a through (|Frank et al.ll2002h 

Mcapt/Mw - i?a/(4a 2 ) = 2 x lO^Wo 2 ,. (2) 

Here a=4.2xl0 12 aiod cm is the orbital separation, aiod=Piod-^30 3 ' Piod is the binary orbital period in units of 
10 days, and M30 is the total mass in units of 30 M Q (we assumed circular orbits). 

• The magnetospheric radius, Rm, at which the pressure of the NS magnetic field (/j, 2 / (8ttR^ s ), with fj, the NS 
magnetic moment) balances the ra m pressure of the inflow ing matter (p w w^). In the case in which RM>R a , the 
magnetospheric radius is given by (|Davies fc Pringlelll981f ) 

R M = 3.3 x 10 10 AC 6 1/6 ^ 1/6 ai^^ 3 cm. (3) 

Here we assumed a non magnetized spherically symmetric wind ([Eisner fc Lamb! Il977f ). with density 5 
Pw(RM)~Mw/(47ra 2 v w ), a dipolar NS magnetic field with (j, 33 =n/10 33 G cm3 > and M-6=^w/10~ 6 M© yr -1 . 
In the following sections we discuss the range of applicability of Eq. (3[ and the regimes in which a different 
prescription for Rm should be used. 

• The corotation radius, R co , at which the NS angular velocity equals the Keplerian angular velocity, i.e. 

R co = 1.7 x 10 10 P S 2 3 /3 cm. (4) 

Here P s3 is the NS spin period in units of 10 3 s. 

Changes in the relative position of these radii result into transitions across different regimes for the NS 
([Illarionov &; Sunvaevi 119751 : IStella etlrtl Il986[ ) . Below we discuss these regimes singularly, and provide a schematic 
representation of each regime in Fig. [TJ Being deter mined primarily by the spin o f the neutron star, the corotation 
radius can change only over evolutionary timescales. Illlarionov &: Sunvaevi (|1975h summarises the different regimes 
experienced by a spinning down NS since its birth, from the initial radio pulsar (i.e. rotation powered) stage, to the 
regime in which mass accretion onto the neutron star surface can take place. On the other hand the accretion radius 
and magnetospheric radius depend on the wind parameters (see Eqs. [T]and[3]), which can vary on a wide range of 

5 We approximated a-Rjy[— a , which is satisfied for a very wide range of parameters. 
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timescales (from hours to months). Therefore, variations in the wind parameters can cause the neutron star to undergo 
transitions across different regimes on comparably short timescales, thus opening the possibility to explain the prop- 
erties of some classes of highly variable X-ray sources through them. In particular, the transition from the accretion 
regime to the propeller regime (and vice versa), across the so-called "centrifugal barrier" , was identified as a likely 
mechanism responsible for the pronounced activity of Be X-ray pulsar transient systems ([Stella et al.l Il986l ). Below 
we summarise the different regimes of a magnetic rotating neutron star, subject to a varying stellar wind, with special 
attention to t he condition under which a "magnetic" (as opposed to "centrifugal" ) barrier inhibits acc retion onto the 
neutron star ([Harding et al.lll992t lRutledgell200U iMori fc Rudermanl 120031 ; IToropina et~aTll2006l I2001D . As it will be 
clear in the following, new motivation for investigating magnetic inhibition of accretion comes from the d iscovery of 
magnetars, neutrons stars with extremely high magnetic fields (^10 14 -10 15 G. lDuncan fc Thompsonlll992l ). 

3.1. Outside the accretion radius: the magnetic inhibition of accretion: Ryi>R a 

We consider here the case in which the magnetospheric radius is larger than the accretion radius 6 . In systems with 
R,M>Ra the mass flow from the companion star in teracts directly with the NS magnetosphe re without significant 
gravitational focusing, forming a bow shock at Rm (jHarding et alj [19921 IToropina et al.|[200U ). A region of shocked 
gas surrounds the NS magnetosphere with density p ps ~4p w and velocity v ps ~v w /4 (the subscript "ps" stands for 
post-shock). These are only rough estimates beca use the shock is very close to the magnetopause and it does not 
satisfy the standard Rankine-Hugoniot conditions (jToropina et al.lfeOOlh . At least in the front part of the shock, i.e. 
in the region around the stagnation point, the whole kinetic energy of the inflowing matter is converted into thermal 
energy, and the expected temperature of the heated gas is T~m p v^,/(3k)~4x 10 7 v| K. Thus, the power released in 
this region is of order 

ishock |i?^/9 w Ww = 4.7 x 10 29 R 2 Mw vla-£ d M^ erg s" 1 (5) 

(Rmi o is the magnetospheric radius in units of 10 10 cm), and is mainly radiated in the X-ray band (jToropina et al.l 
2006). Below we distinguish two different regimes of magnetic inhibition of accretion. 

3.1.1. The super Keplerian magnetic inhibition regime: Ryi>R a , R co 

In this "superKeplerian" magnetic inhibition regime the magnetospheric radius is larger than both the accretion and 
corotation radii (RM>R a , Rco)- Matter that is shocked and halted close to Rm cannot proceed further inward, due to 
the rotational drag of the NS magnetosphere which is locally superKeplerian. Since magnetospheric rotation is also 
supersonic 7 , the interaction between the NS magnetic field and matter at Rm results in rotational energy dissipation 
and thus, NS spin down. In order to derive an upper limit on the contribution of this dissipation to the overall 
luminosity, we assume that the above interaction is anelastic (e.g.. iPerna et a l. 2006), i.e. that matter at Rm is forced 
to corotate. This process releases energy at a rate 

L Bd ~ irR 2 M p w v w {R M n) 2 ~ 3.7 x 10 29 R^M^a-^P^ 2 erg a" 1 , (6) 
which adds to the shock luminosity (Eq. [5J . 

3.1.2. The subKeplerian magnetic inhibition regime: R a <RM<Rco 

If R a <RM<Rco the magnetospheric drag is subKeplerian and matter can penetrate the NS magnetosphere. In this 
"subKeplerian" magnetic inhibition regime, the boundary between the inflowing matter and the magnetosphere is 
subject to the Kelvin- Helm holtz instability (KH I, lHarding et al.iri992T ). The mass inflow rate across Rm resulting from 
the KHI is approximately (|Burnard et al.|jl983l ) 

M K H - 27Ti?MPpsVconv = 27ri?MPps«sh'7KH (Pi/Pe) 1/2 (1 + Pi/ ' Pc)~ l , (7) 

where ?7kh~0.1 is an efficiency factor, v s h is the shear velocity, p\ and p c the density inside and outside the magne- 
tospheric boundary at Rm, respectively. Close to the stagnation point, virtually all the kinetic energy of the wind 
matter is converted into thermal energy (see also § 13.1. 1|) , and the shear velocity is thus dominated by the magne- 
tosphere's rotation, such that v s h=v ro t=27rPg Rm- Away from this region, the tangential component of the wind 
velocity with respect to the NS magnetic field lines increases up to v ps and, in general, the shear velocity and rate at 
which plasma enters the magnetosphere due to the KHI, depends on both v ps and v rot . For the aims of this paper, 
we ado pt v s h=max(v D s, Vrot ), and use mass conservation across the KHI unstable layer to estimate the density ratio 
at Rm (jBurnard et al.l ll983). If matter crossing the unstable layer of height h t is rapidly brought into corotation with 
the NS magnetosphere and free-falls onto the NS, mass conservation implies 

~ R M h t piVs(RM)- (8) 

The height h t of the unstable layer, where matter and m agnetic field coexist, is mostly determined by the largest 
wavelength of the KHI unstable mode (jBurnard et alj [l983). A detailed analysis of this instability is beyond the scope 

6 A similar case was considered also by ILipunovl 1)19921) . "the georotator regime", and byjToropina ct al. (2001), "the magnetic plow 
regime" . 

7 This can be easily seen by comparing v w with ORm, for the value of the parameters used in this section. 
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of this paper. In Eq. [8] we use conservatively h t ~RM, and discuss in Appendix [Al the effect of smaller values of h t . 
Therefore accretion of matter at a rate Mkh onto the NS is expected to release a luminosity of 

Lkh = 3.5 x 10 3 V H i?Mioarod^-6(pi/pc) 1/2 (l + Pi/Pc)' 1 erg s" 1 , (9) 

if v sh =v ps , or 

Lkh = 8.8 x lO^K^p-^R^a-^M^ipi/p,) 1 / 2 ^ + a/aO" 1 erg s~\ (10) 

if v s h=v rot . The values of p\j p c that we use in Eqs. [51 and 1101 are derived numerically from Eq. [H] 

In the subKepleri an magnetic inhibition regim e, plasma penetration inside the NS magnetosphere is sustained also 
by Bohm diffusion (jlkhsanov fc Pustil'nikl[l996[ ). This diffusion, being dependent on the temperature of the plasma 
(rather than velocity), is highest in the region close to the s tagnation point, where the shock slows down the inflowing 
plasma most efficiently (see also § 13. 1 . lj) . In accordance with llkhsanov et al.l (|2001bl ). the maximum inflow rate allowed 
by Bohm diffusion is 

M diff ~ 2irR 2 M p ps V m = 4.5 x lOY^M-eM^^moVod g s" 1 , (11) 

where V m —y/D s/ts is the diffusion velocity, D c ff=(£ckTi(i?M))/(16eB(RM)) the diffusion coefficient, 
Ti(Rm)— m p v w/(3k) the post-shock ion temperature, tg=yj R M / (2GM) the free-fall time, C— 0.1 an efficiency factor 
and m p the proton mass. In the above equation we also approximated the den sity outside Rm, arou nd the stagnation 
point, with p ps (though this might be underestimate by a factor of a few, see lToropina et al"1l2001f K Over the whole 
range of parameters relevant to this work, the diffusion-induced mass accretion rate is orders of magnitude smaller 
than that due to the KHI. Similarly, the contribution to the total luminosity resulting from the shock and anelastic 
drag at the magnetospheric boundary can be neglected in this regime. 

3.2. Inside the accretion radius: Rm<R& 

3.2.1. The supersonic propeller regime: R co <Rm<R& 

Once Rm is inside the accretion radius, matter flowing from the companion star is shocked adiabatically at R a 
and halted at the NS magnetosphere. In the region between R a and Rm this matter redistributes itself into an 
approximately spherical configuration (resembling an "atmosphere"), whose shape and properties are d etermined by 
the in terac tion between matter and NS magnetic field at Rm- This scenario was considered previously bv lDavies et al.l 
(|1979h and lDavies fe Pringld (|1981l ). and we follow here their treatment. These authors demonstrated that hydrostatic 
equilibrium ensues when radiative losses inside R a are negligible (we discuss this approximation in Appendix [B] and 
lCj) and the atmosphere is stationary on dynamical time-scales. Assuming a polytropic law of the form pocp 1+1 /™, the 
pressure and density of this atmosphere are: 

p(R) = Pps v 2 ps [1 + (1/(1 + n))8R a /R] n+1 (12) 

P (R) = p ps [1 + (1/(1 + n))8R a /R] n . (13) 

The value of the polytropic index n depends on the conditions at the inner boundary of the atmosphere, and in 
particular on the rate at which energy is deposited there. 

When the rotational velocity of the NS magnetosphere at Rm is supersonic (see also § 13.1. 1[) , the interaction with 
matter in the atmosphere le ads to dissipation of som e of the star's rotational energy and thus spin-down. In the 
supersonic propeller regime, iDavies fc Pringiel (Jl981) showed that turbulent motions are generated at Rm which 
convect this energy up through the atmosphere, until it is lost at its outer boundary. In this case n—1/2. Accordingly, 
taking into account the structure of the surrounding atmosphere, the magnetospheric radius is given by 

i? M 6 (l + 16i? a /(3i? M )r 3/2 = ^M w a- 2 p- 2 v w . (14) 

This can be approximated by 

Rm 2.3 x W w a%lMI* /9 vt /9 f4i 9 cm. (15) 

Matter that is shocked at ~R a , reaches the magnetospheric boundary at Rm where the interact ion with the NS magnetic 
field draws energy from NS rotation (see also § 13 . 1 . 1[) . According to IDavies fe Pringiel (|1981l ), this contributes 

L sd = 2nR 2 M p{R M )cl{R M ) 5.4 x 10 31 M_ 6 a^Vg ^^(l + 16i? a io/(3i? M io)) 1/2 erg s" 1 (16) 

to the total luminos ity. In the above equation R a i =10 _10 R a and c s (RM)=vg (Rm)=(2GMns/Rm) 1 ^ 2 
(jPavies fc Prin gle 1981). In the supersonic propeller regime the energy released through the shock at R a , 

ishock = ^Rlp w vl ~ 2.6 x 10 29 Rl 10 v 2 8 a w 2 d M_ e erg s" 1 , (17) 



is negligible. 
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3.2.2. The subsonic propeller regime: Rm<R&, Rco, Mw<M\im 

The break down of the supersonic propeller regime occurs when Rm<R CO j i-e., when the magnetosphere rotation 
is no longer supersonic with respect to the surrounding material. The structure of the atmosphere changes and the 
transition to the subsonic propeller regime takes place. Since the rotation of the magnetosphere is sub sonic, the 
atmosphere is roughly adiabatic (ra=3/2), and the magnetospheric radius is given by (|Davies fc Pringlelll981l ): 

i? M 6 (l + 16i? a /(5.RM)r 5/2 = ^M w a- V V- (18) 

This can be approximated by 

R M ~ 2 x 10 10 a% 7 d MZe /7 4 n ' Ihn cm - ( 19 ) 

In the subsonic propeller regime, the centrifugal barrier does not operate because Rm<Rco, but the energy input at 
the base of the atmosphere (due to NS rotational energy dissipation) is still too high for matter to penetrate the 
magnetosphere at a rate M C apt (jDavies fc Prin gle 1981). Nevertheless a fraction of the matter inflow at R a is expected 
to accrete onto the NS, due to the KHI and Bohm diffusion 8 . 

Based on the discussion in § I3.1.2( we estimate the accretion luminosity of this matter by using Eqs. 171 and 1111 (we 
approximate here the surface of interaction between matter and magnetic field with 47rR M ). This gives 

L diff ~ GM m M diB /R m = 4.5 x W^M^a^R^^ 2 ?^ 1 ^ + 16i? alo /(5i? M io)) 3/2 erg s" 1 , (20) 

and 

Lkh — GM^sMkh/ Rns — 

1.8 x 10 35 77KH J P S 3 1 J RMio^-6a 10 d u 8" 1 ( 1 + 16i?aio/(5i?Mio)) 3/2 (,Oi/Pc) 1/2 (l + Pi/XT 1 erg s-\ 

(21) 

for the accretion luminosity arising from matter entering the magnetosphere through Bohm diffusion and KHI, re- 
spectively. For the range of parameters of interest here, Eqs. [5l [20l and [211 show that Lk h dominates. The rota tional 
energy dissipation at Rm (see § 13.2. ip gives a small contribution with respect to Eg. |2"T1 (|Davies fc Pr ingle 1981]): 

L sd = 2irRi lP (R M )n 3 = 2.2 x lO^P^R^M^v^a^l + 16i? a io/(5i? M io)) 3/2 erg s" 1 . (22) 
The subsonic propeller regime applies until the critical accretion rate 

M llm _ 6 = 2.8 x 10 2 P S 3 3 a 2 0dU8 i?^ /2 (l + 16i? alo /(5i?Mio))" 3/2 (23) 

is reached, at which the gas radiative cooling (bremsstralhung) completely damps convective motions inside the at- 
mosphere (see Appendix [C]). If this cooling takes place, direct accretion at a rate M capt onto the NS surface is 
permitted. 

3.2.3. The direct accretion regime: Rm<R&, Rco, M w >Mn m 
If Rm<Rco and matter outside the magnetosphere cools efficiently, accretion onto the NS takes place at the same 
rate M C apt (see Eq. [2]) at which it flows towards the magnetosphere. The corresponding luminosity is 

i acc - GM NS M capt /R NS - 2 x 10 35 M_ 6au >- 4 erg s" 1 ~ 2 x 10 35 M 15 erg s" 1 , (24) 

where Mi5=M capt /10 15 g s _1 . This is the standard accretion regime; the system achieves the highest mass to luminosity 
conversion efficiency. 

4. TRANSITIONS AND PATHS ACROSS DIFFERENT REGIMES 

We explore here the conditions under which transitions across different regimes take place. As emphasised in § |3l 
these transitions occur when the relative positions of Rm, R a , and R co change; we concentrate here on transitions 
that occur in response to variations in the stellar wind parameters. In the following, since Rm depends only weakly 
on the orbital period and the total mass of the system, we fix aiod=l (we explain this choice in §0, and investigate 
variations in the other four parameters: /133, P S 3, vs, and M_6- 

The equations that define the conditions for transitions between different regimes are 9 

R M >R a ^ M-6 < 0.45^ 3W ii a 2 od (25) 
(or equivalently M 15 <0.6^3 3 Vg, see Eqs. [TJ |3]and[2]), for the magnetic barrier; 

Rm > Rco => P s3 < 2.6MI 6 1/ V 1/4 aiod>33 2 - 2.8M^ V V B/ V33 2 ( 26 ) 

8 To our knowledge this is the first application of the KHI to the subsonic propeller regime. 

9 Here we used Eqs. I15l and ll9l for the magnetospheric radius in the supersonic and subsonic propeller regime, respectively. 
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Fig. 2. — Transitions between regimes described in § \3\ for selected values of the parameters 1133, P s 3, and vs (in each panel these 
parameters are indicated in the top right corner). In all cases we fixed aiod=l (see §21). In panels (a) and (b) we fixed /J33 and vg and 
investigated the different regimes in the P s 3-M_g plane. In panels (c), (d), (e), and (f), instead, fixed parameters are P s 3 and /X33, and the 
relevant regimes are shown in the vg-M—e plane. In all cases Eqs. 1251 1261 1271 1281 and IB2I (see Appendix iBt are represente d wi th a solid 
line, a dotted line, a dot-dot-dot dashed line, a dot-dashed line, and a dashed line, respectively. Note that the line from Eq. IB2I is present 
only in panel (e), i.e our treatment of the the supersonic propeller is self-consistent everywhere except for a small region away from the 
range of interest for SFXT sources (see also § [SJ . 
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(see Eq. [3]) if R M > R a , or 

R M > R co => P s3 < 1.8aJ^AC 6 1/3 Us 2/3 ^ 3 - 2M- 1/3 v- 2/3 [i% 3 (27) 

(see Eq. [15]) if RM<R a , for the centrifugal barrier. 
The equation that defines the transition from the subsonic propeller to the direct accretion regime is 

D A r sV-iS/ 21 30/21 60/21 16/21 _ _ 15/21 16/21 ,„ n 

F s3 > 4.5Af_ 6 a 10 ^ v 8 n 33 > ~ 5.5M 15 p^ 7 . (28) 

In Fig. [5] the above equations are represented as lines separating different regimes. In panels (a) and (b) we fixed 
;U33 and vs and investigated the different regimes in the P S 3-M_g plane. In panels (c), (d), (e) and (f), instead, P S 3 
and /i33 were fixed and the relevant regimes shown in the v§-M_6 plane. Below we summarise the different regimes 
that a system attains in response to variations of Af w , in the different panels. 

Panel (a) shows that, for a wind velocity of vg=1.2, a strongly magnetized NS (^33=0.8) can undergo a transition 
between the superKeplerian and subKeplerian magnetic inhibition regimes in response to changes in the mass loss rate 
only for typical spin periods >2000 s. When the mass loss rate reaches M_6^2.4, the direct accretion or subsonic 
propeller regime sets in, depending on whether the spin period is longer or shorter than ~3700 s. For spin periods 
<420 s the direct accretion regime is not attained for the interval of mass loss rates considered in Fig. [2] and only 
transitions between the superKeplerian magnetic inhibition and supersonic propeller regime are expected. 

For lower magnetic fields (/U33=0.01), panel (b) shows that only the supersonic propeller, subsonic propeller, and 
direct accretion regime can be attained. For spin periods in the range ^60-230 s transitions can occur between all these 
three regimes, whereas systems with spin periods longer than ~1300 s and shorter than ~20 s are expected to be in 
the direct accretion and the supersonic propeller regime, respectively. For /Z33=0.01 transitions to the superKeplerian 
and subKeplerian magnetic inhibition regimes cannot take place because the magnetospheric radius is too small to 
exceed the accretion radius. 

In panel (c) (/i33=0.8 and P S 3=1), transitions can occur virtually between all the regimes described in § EI In 
particular, for vs in the range 0.9-1.5 transitions are expected to take place between the superKeplerian magnetic 
inhibition, the supersonic and subsonic propeller, and the direct accretion regime, as the mass loss rate increases 
from M_6~0.1 to M_6~100. For velocities v§<0.9 transitions can occur only between the supersonic propeller, 
the subsonic propeller and the direct accretion regime, whereas transitions to the superKeplerian and subKeplerian 
magnetic inhibition regimes are impeded by the fact that the accretion radius cannot be overtaken by Rm- On the 
contrary, for vg>1.5, the magnetospheric radius is located beyond the accretion radius for any considered value of 
the mass loss rate, and thus transitions can take place only between the superKeplerian and subKeplerian magnetic 
inhibition regimes. Similar considerations apply to panels (d), (e), and (f). 

For a system with /X33=0.01 and P S 3=1 (panel (d)), the superKeplerian magnetic inhibition regime never occurs, 

because the magnetic field is too low and Rm<R co for any 0.1<M_6<100. Instead, the subKeplerian magnetic 
inhibition regime can be attained for high wind velocities (vs=2), because R a ocv~ 2 . 

In panel (e), for ^33=0.8 and P s3 =0.1, the magnetospheric radius is larger than the corotation radius for the entire 
range spanned by M_ 6 . Thus accretion o nto the NS does not ta ke place (note that in the region below the dashed 
line accretion can occur even if R M >.R™. lDavies fc Pringlei ri981) 10 . 

Finally, in panel (f) we show the transitions for a system with /i 33 =0.01 and P s3 =0.1 s. In this case all the regimes 
described in §[3] are present in the figure, similar to the case of panel (c). However, the region corresponding to the 
subsonic propeller regime is larger, such that there is only a modest range of (fixed) velocities for which mass loss rate 
variations in our chosen range (0.1<M_@<100) can cause transitions through all regimes, from the superKeplerian 
magnetic inhibition to the direct accretion regime. As we discuss below, this has important consequences for the 
expected luminosity variations. 

We now compute the luminosity swings for some of the examples discussed a bove in a fash i on similar to what wa s 
done in the context of centrifugally inhibited accretion in NS X-ray transients (|Corbetlll996t ICampana et al.N l998). 
Fig.[3K applies to a system with ^ 33 =0.8, P S 3=1, and v 8 =1.2 (this corresponds to the case P S 3=1 of Fig. [2] panel (a)). 
The lower panel of this figure shows that, for 0.1<M_6<100, the magnetospheric radius crosses both the centrifugal 
(R co ) and magnetic (R a ) barriers. Correspondingly, the system moves from the superKeplerian magnetic inhibition 
regime, to the supersonic and subsonic propeller regime, and, finally, to the direct accretion regime, giving rise to a 
six-decade luminosity swing from ^10 31 to ~10 37 erg s _1 . We note that a large part of this swing (about five decades) 
is attained across the transitions from the superKeplerian magnetic inhibition to the direct accretion regimes, which 
take a mere factor of ~5 variation of M w . 

In the presence of a standard NS magnetic field (10 12 G), Fig. [3)3 shows that such abrupt luminosity jumps are not 
expected for a very slowly rotating (1000 s) NS (the other system parameters are the same as those of Fig. [3K), since 
the magnetospheric radius is smaller than both R a and R co , for any reasonable value of M w . Therefore, the direct 
accretion regime applies, with the the luminosity proportional to M w . 

10 This is because, in this region, the mass flow rate is so high that Eq, IB2l is violated (the supersonic propeller is no longer self-consistent), 
and convective motions are damped by radiative cooling. 



Magnetars and SFXTs 



9 




Fig. 3. — (A) Upper panel: Variation of the luminosity through different regimes, as a function of the mass loss rate from the companion 
star. In this case the parameters of the model are fixed at /i33=0.8, P s 3=l, and vs=1.2. Lower panel : Relative position of the magnetospheric 
radius, Rm (solid line), with respect to the accretion radius R a (dotted line), and the corotation radius R co (dashed line), as a function of 
the mass loss rate from the companion star. 

(B) Upper panel: Variation of the luminosity through different regimes, as a function of the mass loss rate from the companion star. In 
this case the parameters of the model are fixed at /i33=0.001, P s 3=l, and vs=1.2. Lower panel: Relative position of the magnetospheric 
radius, Rm (solid line), with respect to the accretion radius R a (dotted line), and corotation radius R co (dashed line), as a function of the 
mass loss rate from the companion star. 



In Fig. 2] we show the transitions for a system with /i33=0.01 and P S 3=0.1. The wind velocity is vg=1.2 in Fig.[3J\., 
and V8=2.2 in Fig. [3J3 (see also panel (f) of Fig. [5]). These two figures show that, for sub-magnetar fields, a 100 s 
spinning NS can undergo a transition across the magnetic barrier (besides the centrifugal barrier), for suitable 
parameters (a high wind velocity in the case at hand). Such transitions take place over a more extended interval of 
mass loss rates. For instance Fig. @}3 shows that an increase by a factor ~100 in the mass loss rate is required, in this 
case, to achieve a factor ~10 5 luminosity swing comparable with the magnetar case of Fig. 

In order to illustrate further the role of the magnetic field, spin period and wind velocity, we show in Fig. [5] the way 
in which the transitions across regimes take place, by holding two of the above variables fixed and stepping the third 
variable. Figure [SJ^ shows the effect of increasing the value of vs from 1 to 1.8 (in turn resulting in a decrease of the 
accretion radius), in a system with ^t 33 =0.8 and P S 3=1. 

The behaviour of the luminosity changes mainly because a different set of regimes is involved in each case. For Vg=l 
(solid line) the system passes through the superKeplerian magnetic inhibition regime, the supersonic, and subsonic 
propeller regime, finally reaching the direct accretion regime at M_6— 6. In the case vg=1.2 (dotted line), the transition 
to the supersonic propeller shifts towards higher mass loss rates, such that superKeplerian magnetic inhibition applies 
up to M_6<2. Further increasing the wind velocity to vg=1.4 (dashed line), the system first undergoes a transition to 
the subsonic propeller at Af-6— 10, bypassing the supersonic propeller (this is because for M_6— 10 the magnetospheric 
radius defined by Eq. [TJ]is smaller than the corotation radius). As the mass loss rate increases further, the direct 
accretion regime sets in for M_6— 20. In the case vg=1.8 (dot-dashed line), the corotation radius exceeds the accretion 
radius. The system is thus in the superKeplerian magnetic inhibition regime for Af_g<20, while for higher mass loss 
rates the subKeplerian magnetic inhibition regime applies and the luminosity is dominated by accretion through the 
KHI. A transition to the direct accretion regime is expected for mass loss rates M_6>100. 

Figure [5)3 shows the transitions across different regimes for selected values of the magnetic field, in a system with 
P S 3=1 and vg=1.2. For the lowest magnetic field in this figure (/i33=0.01, solid line), the magnetospheric radius is 
smaller than both the accretion and corotation radius for the whole range of mass loss rates spanned in the figure, and 
the system is always in the direct accretion regime. In the case //33=0.1 (dotted line), the system is in the subsonic 
propeller for M-q<1, while for higher mass inflow rates the direct accretion regime applies. For /i33=0.5 (dashed line) 
the luminosity behaviour becomes more complex, as the system goes through the superKeplerian magnetic inhibition 
regime, the supersonic and subsonic propeller regime, and eventually reaches the direct accretion regime at M__6— 8. 
For //33=0.8 the same sequence of transitions applies, with the entire luminosity swing taking place over a smaller 
interval of mass loss rates. 

Finally, in Fig. [5p we show the effects of increasing the spin period in a system with /Z33=0.1 and vs=1.2. For the 
lowest spin period considered here (P S 3=0.01, solid line), the system remains in the supersonic propeller regime for 
the whole range spanned by M_e- In the case P S 3=0.5 (dotted line) transitions occur from the supersonic propeller 
regime (M_6<0.7), to the subsonic propeller (0.7<M_6<4) and then to the direct accretion regime (M_6>4). By 
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Fig. 4. — (A) Upper panel: Variation of the luminosity through different regimes, as a function of the mass loss rate. In this case the 
parameters of the model are fixed at /i33=0.01, P s 3=0.1, and vs=1.2. Lower panel: Relative position of the magnetospheric radius, Rm 
(solid line), with respect to the accretion radius R a (dotted line), and corotation radius R co (dashed line), as a function of the mass loss 
rate from the companion star. 

(B) Same as (A) but for /i 33 =0.01, P s3 =0.1, and v s =2.2. 



further increasing P S 3 to a value of 1 (dashed line), the system goes through the subsonic propeller and the direct 
accretion regime, while the supersonic propeller regime does not occur due to the longer spin period as compared to 
the previous case. For P S 3=6 (dot-dashed line) the spin period is so high that the direct accretion regime applies for 
any reasonable value of Af_g. 

The above results show that over a range of values of the key parameters /133, P S 3, and vs, large luminosity swings can 
be achieved with comparatively modest changes in the mass loss rate, as the neutron star undergoes transitions from 
one regime to another. More generally, these transitions result from changes in the relative position of the accretion, 
magnetospheric, and corotation radii, reflecting short term variations of the wind velocity and mass loss rate, the only 
parameters that can vary on shorter timescales than secular. Therefore, transitions between different regimes take 
place once the source parameters are such that the NS straddles the centrifugal barrier (i.e. Rm— Rco, when R a >R co ) 
or the magnetic barrier (i.e. Rm— R a , when R a <R co ) or both (i.e. Rm— Rco— Ra)- The centrifugal barrier applies to 
relatively short spin periods. It is well known that the longer the spin period of a transient neutron star, the higher 
its magnetic field must be for the centrifugal barrier to operate (see Eqs. l26l and l27|). In particular, for periods of 
hundreds seconds, or longer, magnetic field strengths of >10 13 -10 14 G are required 11 . On the other hand, we have 
shown that neutron stars with even longer spin periods and magnetar-like fields are expected to undergo transitions 
across the magnetic barrier and thus are expected to have an inherently different "switch off" mechanism than short 
spin period systems. A necessary condition for this is that Rc >R a , which translates into Ps3>3v^~ 3 . The magnetic 

and centrifugal barriers set in (nearly) simultaneously (i.e. R a ~R co ~R M ) for fi 33 ^0.3P 7 s ^ 6 M^ 2 . 
Taking into account of all the examples discussed in this section, we conclude that: 

1. Long spin period systems (P S 3<^1) require magnetar-like B-fields (a*33^0.1) in order for a large luminosity swing 
(~10 5 ) to arise from modest variations in the wind parameters (e.g. a factor ~5 in M-e). These luminosity 
swings might result from transition across different regimes through both the centrifugal and magnetic barriers. 

2. Shorter spin period systems (P S 3<?Cl) must posses lower magnetic fields (/X33<C0.1) for similar transitions to take 
place. Somewhat larger variations in the wind parameters are required in order to achieve similar luminosity 
swings to those of the long period case, and transitions between different regimes occur in most cases through 
the centrifugal barrier. 

3. Few or no transitions are expected for systems with either high magnetic fields and short spin periods, or systems 
with lower magnetic fields and long spin periods. In the first case the centrifugal barrier halts the inflowing 
matter at Rm and accretion does not take place; such systems might thus be observable only at very low (X- 
ray) luminosity levels (~10 32 -10 33 erg s _1 ). In the second case Rm<R co for a wide range of wind parameters, 
accretion can take place, and a high persistent luminosity is released (~10 35 -10 37 erg s _1 ). 



5. APPLICATION TO SFXT SOURCES 

In this section we propose that transitions across different regimes caused by relatively mild variations of the wind 
parameters are responsible for the outbursts of SFXTs. In consideration of the wide range of spin periods inferred for 



11 This point was already noted in Stella ct al. (1986). 




Fig. 5. — Variations of the luminosity through different regimes, as a function of the mass loss rate, for different sets of parameters. 

(A) : In this case we fixed ^33=0.8 and P s 3=l. The solid line shows luminosity variations for vg=l, the dotted line is for vg=1.2, and the 
cases V8=1.4 and V8=1.8 are represented with a dashed and dot-dashed line, respectively. 

(B) : Here the fixed parameters are P s 3=l and vg=1.2. The different curves correspond to ^33=0. 01 (solid line), 0.1 (dotted line), 0.5 
(dashed line), and 0.8 (dot-dashed line). 

(C) : In this case we fixed ^33=0. 1 and vs=1.2. Luminosity variations across different regimes are shown for P s 3=0.01 (solid line), 0.5 
(dotted line), 1 (dashed line), and 6 (dot-dashed line). 



SFTXs, the outbursts of these sources are expected to result from the openi ng of the magnetic barrier in very long 
spin period systems and the centrifugal barrier in all other systems (see also iGrebenev fc Sunvaevll2007t ISidoli et al.1 
|2007|) . More crucially, we conclude that slowly spinning SFXTs should host magnetars, independent of which of the 

two mechanisms is responsible for their outbursts; 

As a case study we consider IGR J17544-2619 (|Sunvaev et al.l[20 03i) , a SFXT observed by Chandra during a complex 
transition to and from a ~1 hour-long outburst, yielding t he first de t ailed characterization of a SFXT light curve. 
The spin period of IGR J17544-2619 is presently unknown. Iln't Zand! (|2005l ) showed that four different stages, with 
very different luminosity levels, could be singled out during the Chandra observation: (a) a quiescent state with 
Lx— 2xl0 32 erg s _1 , (b) a rise stage with Lx— 1.5xl0 34 erg s _1 , (c) the outburst peak with Lx— 4xl0 37 erg s _1 , and 
(d) a post-outburst st age (or "tail" ) with L x— 2xl0 36 erg s _1 (see panel (a) of Fig.[6l these luminosities are for a source 
distance of ^3.6 kpc. lRahoui et aLl l2008). The maximum luminosity swing observed across these stages was a factor 
of>6.5xl0 4 . 

Motivated by the evidence for >1000 s periodicities in XTEJ1739-302 and IGR J16479-4514 (see §E), we discuss 
first the possibility that IGR J17544-2619 contains a very slowly spinning neutron star. We use ^33=1, P S 3=L3, 
vg=1.4, and show in Fig. [U[b) the different regimes experienced by such a neutron star as a function of the mass 
loss rate. For M_ 6 <20 the above values values give Rm>R & and Rm>R co , such that superKeplerian magnetic inhi- 
bition of accretion applies. The expected luminosity in this regime, ^10 31 erg s _1 , is likely outshined by the X-ray 
luminosity of t he supergiant star (the companion star's lum inosity is not shown in Fig. [6j but it is typically of order 
~I0 32 erg s -1 . ICassinelli et al.llf 981t iBerghoefer et aLll2001h . We conclude that the lowest emission state (quiescence) 
of IGR J17544-2 619 can be explained in this way, with the companion star dominating the high energy luminosity 
(jln't Zand! 120(151 ). The rise stage is in good agreement with the subKeplerian magnetic inhibition regime, where the 
luminosity (~T0 34 erg s _1 ) is dominated by accretion of matter onto the NS due to the KHI. The uncertainty in the 
value of h translates into an upper limit on the luminosity in this regime which is a factor of ~10 higher than that 
given above (see § 3 and Appendix |A|). 

During the outburst peak the direct accretion regime must apply at a mass loss rate of Af_g=500. In this interpreta- 
tion direct accretion must also be at work in the outburst tail at M_6^3, where a slight decrease in M w would cause the 
magnetic barrier to close and the source to return to quiescence. According to this interpretation, if IGR J 17544-26 19 
has a spin period of >I000 s, then it must host a magnetar. 

Panel (c) of Fig. [6] shows an alternative interpretation of the IGR JI7544-26I9 light curve, where we fixed /i33 =0.08, 
P S 3=0.4, and Vg=l. For this somewhat faster spin (and lower magnetic field), the luminosity variation is mainly driven 
by a transition across the centrifugal barrier (as opposed to the magnetic barrier). In this case, the quiescent state 
corresponds to the supersonic propeller regime (M-e <0.6), the rise stage to the subsonic propeller (0.6<M_6<2), 
while both the peak of the outburst and the tail take place in the direct accretion regime at M_e=200 and M_g=I0, 
respectively. 

Assuming an even faster NS spin period for IGR J17544-2619, a weaker magnetic field would be required. In panel 
(d) of Fig. [6l we show the results obtained by adopting ^33 =0.001, P S 3=0.01, and v§=2. The ~I0 34 erg s _1 luminosity 
in the subsonic propeller regime compares well with the luminosity in the rise stage. However, the luminosity of 
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the supersonic propeller regime is now significantly higher than the quiescence luminosity of ~10 32 erg s 1 (this is 
consequence of the higher value of M w for which the supersonic propeller regime is attained in this interpretation). 
We note that, the whole luminosity swing takes place for a wider range of mass loss rates, and the outburst peak 
luminosity requires Af_ 6 ~3000, an extremely high values even for an OB supergiant. 

Interpreting the properties of IGR J17544-2619 in terms of a NS with a spin periods <C100 s is more difficult. For 
instance, for the subsonic propeller regime to set in, the mass loss rate corresponding to the transition across Rm — B co 
must be lower than the limit fixed by Eq. 1231 If instead the transition takes place at higher mass loss rate, the system 
goes directly from the supersonic propeller to the direct accretion regime (or vice versa), bypassing the subsonic 
propeller: therefore, the rise stage would remain unexplained. Since fast rotating NSs require lower magnetic fields for 
direct accretion to take place while in outburst, Eq. [23] is satisfied only for very high wind velocities (vg>2-3). On the 
other hand, an increase by a factor of ~2 in the wind velocity (with respect to the longer spin period solutions) would 
give a substantially lower M cap t, such that the subsonic and the direct accretion regime luminosities fall shortwards 
of the observed values (unless unrealistically high mass loss rate are considered). 

Based on the above discussion, we conclude that IGR J17544-2619 likely hosts a slowly rotating NS, with spin period 
>100 s. Whether the magnetic barrier or the centrifugal barrier sets in, causing inhibition of accretion away from the 
outbursts, will depend on whether the spin period is longer or shorter then ~1000 s. We note that IGR J16418-4532, 
a ~1240 s pulsating source with a 3.7d orbital period, displayed short duration flares similar to those of SFXTs and 
thus might be considered a candidate for hosting an accreting neutron star with a magnetar-like field. However there 
is no clear evidence yet that I GR J16418-4532 is a t ransient source, since the very low st ate revealed with SWIFT 
might well be due to an eclipse (|Tomsick et alj|2004 IWalter et al.|[200l ICorbet et al.ll2006f ). 

As another example we discuss the case of IGR J16465-4507, a SFXTs with a spin period of 228 s. The luminosity 
behaviour of this sourc e is still poorly known. An outburst at 5x 10 36 erg s _1 was observed with INTEGRAL (assuming 
a distance of 12.5 kpc. iLutovinov et al.1 12005: Smith 2004), which did not detect the source before the outburst down 
to a level of 5xl0 35 erg s -1 . About a week later, XMM-Newton revealed the source at 5xl0 34 erg s _1 and discovered 
the 228 s pulsations (jLutovinov et alJfeOOSt IZurita Heras fc Walterl l2004) . If the direct accretion regime applied all the 
way to the lowest luminosity level observed so far, then an upper limit of /Z33 = 0.004 would be obtained by imposing 
that the neutron star did not enter the subsonic propeller regime (see Eq. I28[) . On the other hand, if the luminosity 
measured by XMM-Newton signalled that the source entered the subsonic propeller regime, while direct accretion 
occurred only during the outburst detected by INTEGRAL, then a considerably higher magnetic field of fi^=0.07 
would be required. 

The above discussion emphasizes the importance of determining, through extended high sensitivity observations, 
the luminosity at which transitions between different source states occur, in particular the lowest luminosity level for 
which direct accretion is still at work. In combination with the neutron star spin, this can be used to infer the neutron 
star magnetic field. Alternatively accretion might take place unimpe ded at all luminosity levels of SFXTs, a possibilit y 
which requires a very clumpy wind as envisaged in other scenarios (jNeeueruela et alj 2008; Zu rita Heras et alJl2007j ). 
In this case the neutron star magnetic field can be considerably lower than discussed here. More extensive studies of 
these sources (and, by analogy, other SFXTs) are clearly required. 

6. DISCUSSION 

If the centrifugal barrier operates in IGR J17544-2619, then the activity of this source (and by extension that of 
SFXTs with simila r properties) should parallel that of long spin period X-ray pulsar transients with Be star companions 
(jStella et al.lll986l ). One crucial difference between the two classes, namely the duration of the outbursts, might well 
result from the presence in Be systems of an accretion disk mediating the flow of matter outside the NS magnetosphere. 
In fact, while in Be systems the star's slow equatorial wind has enough angular momentum to form such a disk, the 
supergiant 's wind in SFXTs is fast and possesses only little angular momentum relative to the NS. In the absence of 
a disk, variations in the wind parameters take effect on dynamical timescales, whereas in the presence of a disk they 
are smoothed out over viscous timescales. 

If the spin period is sufficiently long in IGR J17544-2619 and other SFXTs, the onset of the magnetic barrier will 
inhibit accretion. While steady magnetic inhibition of accretion is familiar, e.g. from the earth magnetosphere - 
solar wind interaction, transitions in and out of this regime have not yet been observed, to the best of our knowledge. 
Therefore, very long spin period SFXTs might provide the first opportunity to study transitions across such a magnetic 
barrier. Irrespective of whether the centrifugal or magnetic barrier operates in IGR J17544-2619, a long spin period 
would imply a high magnetic field, comparable to those inferred for magnetar candidates. 

Sce narios involving magnetars with spin perio ds well above 1000 s have been considered in several studies ([Rutledgd 
120011: iMori fc Rudermani 120031 : iToropina et ail 120011: iLiu fc Yanl l2006t IZhang et alj I2004D . Moreover a few known 
accr eting X-ray pul sars with unusually long spin periods have already been proposed as magnetar candidates (see 
e.g.. Ilkhsanovl [20071 and references therein). However these sources display different properties from SFXTs: some 
are persistent sources; others display week to month-long o utbursts; the high spin up measured in IGR J16358-4726 
testifies that the accretion flow is likely mediated by a disk (|Patel et al.ll2007l ). Therefore some of the features of the 
model discussed in this paper would not be applicable t o these sources. 

In the context of wind- fed HMXBs, Zha ng et al.l (|2004l ) pointed out that magnetars might be hosted in binary systems 
with relatively short orbits (~1-100 d) and long pulse period (>10 3 -10 5 s). By using evolutionary calculations, these 
authors showed that magnetars would be easily spun-down to such long spin periods by the interaction with the wind 
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Fig. 6. — Application of the model described in §[3]to the IGR J17544-2619 transition from quiescence to outburst. 
Panel (a): A schematic representation o f the Cha ndra light curve of IGR J17544-2619 obtained by using the segments that were not affecte d 
by pile-up (see both panels of figure 2 in In't Zand 2005). The different luminosity stages are clearly visible. According to In't Zand (2005), 
the count rates on the y-axis correspond to 2 X 10 32 erg s _1 , 1.5 X 10 34 erg s — x , 2 X 10 36 erg s — 1 , and 4 X 10 37 erg s _1 , in the quiescence 
state, the rise state, the tail, and the peak of the outburst, respectively. 

Panel (b): Interpretation of the quiescence to outburst transition of IGR J17544-2619 in terms of the magnetic barrier model. The dotted 
horizontal lines mark the luminosity that divide the different regimes. The parameters of the model are fixed at P a 3=1.3, fg=1.4, and 
M33=l- 

Panel (c): Interpretation of the quiescence to outburst transition of IGR J17544-2619 based on the centrifugal barrier model. The parameters 
of the model are fixed at ^33=0. 08, P s 3=0.4, and 1)8=1. 

Panel (d): Same as panel (c) but here the parameters of the model are fixed at £133=0. 001, P s 3=0.01, and i)g=2. 



of the companion star, in less than 10 6 yr. The systems we considered in this paper would likely result from a similar 
evolutionary path. 

Once a system approaches the spin period that is required for direct accretion to occur, the short timescale (~ hours) 
erratic variations that characterize a supergiant's wind will cause transitions across different NS regimes. Insofar as 
the average wind properties evolve only secularly, the NS spin will then remain locked around such period, alternating 
spin-up intervals during accretion and spin-down intervals during quiescence, when a ccretion is inhibite d. 

The relevant spin-up timescale during accretion intervals is approximately (see e.g. lFrank et al.| [2002) 

T SU = -Pspin/Pspin = M / '(M capt i) 8 X 10 2 / 45 «t Aod^Mfy 1 (29) 

where l=27rR 2 /(4P or b) the specific angular momentum of wind matter at the accretion radius, I=10 45 ±45 g cm 2 the 
NS moment of inertia, and Mi7=M capt /10 17 g s" 1 (this corresponds to an outburst luminosity of ^10 37 erg s _1 ). 

A rough estimate of the spin-down timescale is obtained by assuming that most of the quiescence luminosity draws 
from the rotational energy of the NS {D avies k Pringldfl98lh : this gives 

T sd = JW-Pspin = m 2 /L X = 1.3 X 10 2 / 4 5i 3 l lp s3 2 (30) 

where L31 is the luminosity produced by the interaction between the magnetosphere and the wind in units of 
10 31 erg s _1 . It is apparent that, for magnetar-like fields and long spin periods the above timescales are much 
shorter than the lifetime of the supergiant's strong wind phase. Therefore secular changes in the wind parameters 
and/or the NS magnetic field will be easily tracked by the NS spin. Moreover since spin-up and spin-down take place on 
comparable timescales, it is to be expected that spin-up intervals during outbursts are compensated for by spin-down 
during quiescent intervals of comparable duration. In other words accretion state of long spin period SFXTs would 
be expected to have a duty cycle of order ^0.5 or higher. As we discussed in § [2j evidence for high values of the duty 
cycle is gradually emerging from high sensitivity observations of SFXTs. 

On the other hand spin-down may also occur during the accretion intervals as a results of velocity and density 
gradients in the supergiant's wind that lead to temporary reversals of the angular momentum of the captured wind 
relative to the neutron star (note that persistent w ind accreting X -ray pulsars in HMXBs have long been known to 
alternate spin-up and spin-down intervals, see e.g. iHenrichsl [1998) . This would tend to favor spin-down, such that 
the NS spin might gradually evolve longwards of the spin period that is required to inhibit direct accretion and a 
transient source becomes a persistent one. Interestingly, this might apply to 2S 0114+65, a persistent X-ray pulsar 
with a luminosity of ~10 36 erg s _1 , which displayed variations by a factor of <10. Since its spin period is extremely 
long (~ 10 4 ) i?M is smaller than the corotation radius and the NS accretes continuously from the relatively weak 
supergiant companion's wind (see also lLi fc van den HeuveHH999T l. 

The SFXTs population might thus represent those supergiant HMXBs systems that have not (yet) evolved away 
from the spin period at which transitions across different NS state can take place. 
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We note that for SFXTs to host magnetars their dipole magnetic field must retain values in the At33=0. 1-1 range 
for a few 10 6 yr, i.e. the typical timescale from the formation of the NS to the onset of the supergiant's strong 
wind. Presently known magnetar candidates (soft gamma re peaters, SGR, and anomalo us X-ray pulsars, AXP) have 
estimated ages in the 10 4 -10 5 yr range (for a review see e.g. IWoods fc T hompson 2006f). Little is known on the long 
term evolution of the magnetic field of magnetars and different models have been proposed which lead to different 
predictions. We note that if the irrot ational mode of ambip olar diffusion dominates the B-field decay, /Lt33=0.1-0.3 can 
be expected for ages of a few 10 6 yr (|Hevl fc Kulkarnilll998h . 

According to our proposed scenario the combination of a long spin period and a very large luminosity swing is 
indicative of the presence of a magnetar. This can be f urther corroborat ed through other magnetars signatures, such 
as e.g. proton cyclotron features in th e X-ray spectrum (|Zane et al.ll2~001l) or sporadic subsecond bursting activity such 
as that observed in AXPs and SGRs (|Gavriil et al.ll2002ft . 

Finally we remark on the orbital period of SFXTs: in all regimes described in § the luminosity scales with 
LxocP or j^ 3 v~ 4 . For our fiducial wind parameters, orbital periods of tens of days are required for the transition 
between lo w and high luminosity state s to occur at ~10 36 erg s _1 , a typical luminosity for the onset of SFXT 
outbursts (jWalter fc Zurita Herasl [2007h . This is why t hroughout this paper w e scaled our equations by P or b=10 d 
and used the same value in the examples of Figs. iNegueruela et al.l (|2008l ) showed that orbital periods around 
^10 d are compatible with the NSs being embedded in the clumpy wind from the supergiant companions, rather than 
in a quasi-continuous wind. In this case SFXT outburst durations might be associated with the transit time of a clump 

r out ~ a/v w = 4.2 x 10 i a wd v^ 1 s, (31) 

in reasonable agreement with the observed durations of individual flares (see §[2]). 

We conclude that dumpiness of the stellar wind, an often-used concept for interpreting SFXT activity, applies to 
the gating scenarios described here as well. The main advantage of introducing a gating mechanism rests with the 
possibility to model the very large luminosity swings of SFXTs with much milder density (or velocity) contrasts in the 
wind. 

7. CONCLUSIONS 

In this paper we reviewed the theory of wind accretion in HMXBs hosting a magnetic neutron star with a supergiant 
companion, and considered in some detail the interaction processes between the inflowing plasma and the magne- 
tosphere that are expected to take place when direct accretion onto the neutron star surface is inhibited. We then 
applied this theory to SFXTs and showed that their large luminosity swings between quiescence and outburst (up 
to a factor of ^10 5 ) can be attained in response to relatively modest variations of the wind parameters, provided 
the system undergoes transitions across different regimes. Expanding on earlier work, we found that such transitions 
can be driven mainly by either: (a) a centrifugal barrier mechanism, which halts direct accretion when the neutron 
star rotation becomes superKeplerian at the magnetospheric radius, a mechanism that has been discussed extensively 
in Be star X-ray transient pulsars, or (b) a magnetic barrier mechanism, when the magnetosphere extends beyond 
the accretion radius. Which mechanism and wind interaction regime apply will depend sensitively on the NS spin 
period and magnetic field, besides the velocity and mass loss rate in the supergiant's wind. In particular, the mag- 
netic barrier mechanism requires long spin periods (>1000 s) coupled with magnetar-like fields (>10 14 G). On the 
other hand, magnetar-like fields would also be required if the centrifugal barrier set in at relatively high luminosities 
(> 10 36 erg s _1 ) in neutron stars with spin periods of hundreds seconds. 

Evidence has been found that the spin periods of a few SFXTs might be as long as 1000-2000 s. Motivated by this, 
we presented an interpretation of the activity of IGR J17544-2619 (whose spin period is unknown) in terms of the 
magnetic barrier by a 1300 s spinning neutron star and showed that the luminosity stages singled out in a Chandra 
observation of this source are well matched by the different regimes of wind-magnetosphere interaction expected in 
this case. We discussed also an interpretation of this source based on the centrifugal barrier and a slightly shorter 
spin period (400 s), which reproduced the luminosity stages comparably well. We emphasise that in both solutions 
the required magnetic field strength (>10 15 G and >8xl0 13 G, respectively) are in the magnetar range. 

While the possibility t hat magnetars are hosted in binary sy stem with supergiant companions has been investigated 
by several authors (e.g., IZhang et al.l l2004t iLiu fc Yanl 120061 ) . clear observational evidence for such extremely high 
magnetic field neutron star in binary systems is still missing. According to the present study, long spin period SFXTs 
might provide a new prospective for detecting and studying magnetars in binary systems. 

EB thanks the CEA Saclay, DSM/DAPNIA/Service d'Astrophysique for hospitality during part of this work. MF 
acknowledges the French Space Agency (CNES) for financial support. We would like to thank the referee for useful 
comments. This work was partially supported through ASI and MUR grants. 

APPENDIX 

ON THE HEIGHT OF THE KHI UNSTABLE LAYER 

In this section we expand on the approximation Ii^Rm, introduced previously in § 13.1.21 As discussed by 
iBurnard et al.l (|1983l ). the height h t of the layer where matter and magnetic field coexist due to the KHI, is mostly 
determined by the largest wavelength unstable mode of the KHI itself. These authors suggested that 

h t s hH s , (Al) 
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where h is a factor of order ~1 5 and H s =R^j kb T(Rm)/(GM m p ) is the scale height of the magnetosheath (note 
that H s is roughly of the same order as Rm— 10 10 cm for T~10 8 K). In the subsonic propeller regime, Eq. IAH might 
be a reasonable assumption due to the presence of an extended atmosphere around the magnctosphcric boundary; 
on the contrary it cannot be used in the subKeplerian magnetic inhibition regime, where matter flowing toward the 
NS is shocked very close to the magnetospheric boundary. Despite all these uncertainties, we show below that the 
assumption h t =h Rm, with h~l, gives a conservative estimate of the mass accretion rate due to the KHI. Using Eq. [8] 
and the definition v conv =r/KHV s h(/Oi//o c ) 1 ^ 2 (l+Pi/Po)^ 1 , we derive the equations that define the ratio of the densities 
inside (p;) and outside (p e ) the magnetosphere. These are 

( Pi / Pc ) 1/2 (l + Pi/Pe) = O^KH^^o^ 1 , (A2) 

if v sh =v rot , and 

(Pi/pc) 1 / 2 ^ + pi/pe) = O.It^h/V 1 ^ 2 ,^, (A3) 

if v s h=v ps . Equations I A2I and l A3I show that Pi/p e is an increasing function of h _1 (for fixed values of vs, P S 3 and Rmio)- 
Therefore, being v conv oc(pi/p c ) 1 / 2 (l+p i /p c ) _1 and AiKH°cv conv , the KHI rate of plasma entry inside the magnetosphere 
is also an increasing function of h" 1 (provided that Pi<p c )- Since the KHI unstable layer is located inside the NS 
magnetosphere, the maximum height attainable is h t =RM, an d thus the approximation used in § 13.1.21 and § 13.2.21 
gives a lower limit on the mass flow rate controlled by the KHI in both the subKeplerian magnetic inhibition and the 
subsonic propeller regime. 

Note also that the above lower limit does not violate the st ability condition of the quasi-static atmosphere in the 
subsonic propeller regime. In fact, following Ikhsanovj (|2001d ) . this atmosphere can be considered quasi-static if the 
relaxation time scale of the envelope is less than the drift time scale of the mass flow crossing the magnetospheric 
boundary. In our case matter penetration inside the NS magnetosphere in the subsonic propeller regime is mostly 
provided by the KHI, and the above condition translates into A/KH^(RM/Rco) 3 ^ 2 Afcapt, which is satisfied for a wide 
range of parameters (see Figs. El IH and [5]). 

An upper limit to the KHI mas s flow rate can b e obt ained by assuming pi/p =l (a solution adopted, for example, 
bv lRutledgell2001fh According to iBurnard et all (jl983l ). the density Pi can be increased until the thermal pressure 
inside the magnetosphere piocpi is comparable to the magnetic pressure p m =B 2 (RM)/(87r). When this limit is reached, 
an instability occurs on the lower surface of the unstable layer that increases h t until pi<p m is restored. Using the 
post-shocked gas temperature at Rm (see § E|), it is shown that pi/p e ~l does not violate the condition p;<p m (at 
least in the subKeplerian magnetic inhibition regime). Therefore, even though we restricted ourself to the lower limit 
h=l, the upper limit on Mkh might be attainable in some instances. The KHI would then provide matter penetration 
inside the magnetosphere at a rate ~M capt , thus allowing almost all the captured matter to accrete onto the NS. A 
detailed calculation of KHI accretion in the subKeplerian magnetic inhibition and subsonic propeller regimes will be 
reported elsewhere. 

RADIATIVE LOSSES IN THE SUPERSONIC PROPELLER 

As shown bv lDavies fc Pringli (|198lD . the treatment of the supersonic propeller regime (see § 13.2.10 is self consistent 
only if the energy input at the base of the atmosphere, due to the supersonic rotating NS magnetosphere, is larger 
than radiative losses within the atmosphere itself. The range of validity of this assumption can be determined by using 
the convective efficiency parameter 

r = m%v t (R)t bI (R)R-\ (bi) 

where 9Jtt=v t (R)/c s (R) is the Mach number, v t and c s are the turbulent and sound velocity, and tb r =2xl0 11 
T 1 / 2 m p p _1 (R) s is the bremsstralhung cooling time. For most of the NS rotational energy dissipated at the mag- 
netospheric radius to be convected away through the atmosphere's outer boundary, T should be >1 across the entire 
envelope. Taking into account that in the supersonic propeller regime vt— c s , c s ^Vff, T~Tg and p is given by Eq. 1131 
one finds TocR -3 ' 2 . Therefore the said requirement is satisfied when r(R a )>l, i.e. 

M_ 6 < 2.2 x 10 2 ^a 2 0d . (B2) 

For mass loss rates larger than the above limit, radiative losses are not negligible and the trea tme nt used in § 13.2.11 
for the supersonic propeller regime is no longer self-consistent. We checked that the l i mit of Eq. IB2I is never exceeded 
in the cases of interest. We note that a similar value was also derived by llkhsanovl (|2002| ). but his limit is a factor 
10 larger than ours. This might be due to the fact that llkhsanovl (|2002| ) used p w instead of p ps (l+16/3) 1 / 2 in the 
expression for the matter density at R a . 

RADIATIVE LOSSES IN THE SUBSONIC PROPELLER 

A similar calculation to that in Appendix IB"! in the subsonic propeller regime shows that TocR 1 / 2 , and radiative losses 
are negligible if F(Rm)>1, i.e. 

Af_ 6 < 2.8 x 10 2 P S 3 3 a 2 odU8 i4 /2 (l + 16iW(5i4iio)r 3/2 , (CI) 
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that is the same value as that in Eq. [23] (a somewhat different value was obtained by llkhsanovl (|2001af ) assuming a 
density ~p w at Rm instead of /9 ps (l+16RM/(5R a )) 3 / 2 )- 
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